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Abstract. We describe a pseudo-3D photoionization code, NEBU_3D and its associated visualization tool, VIS_NEB3D, 
which are able to easily and rapidly treat a wide variety of nebular geometries, by combining models obtained with a ID 
photoionization code. We also present a tool, VELNEB_3D, which can be applied to the results of ID or 3D photoionization 
codes to generate emission line profiles, position-velocity maps and 3D maps in any emission line by assuming an arbitrary 
velocity field. As examples of the capabilities of these new tools, we consider three very different theoretical cases. The first 
one is a blister HII region, for which we have also constructed a spherical model (the spherical impostor) which has exactly 
the same Hj8 surface brightness distribution as the blister model and the same ionizing star. The second example shows how 
complex line profiles can be obtained even with a simple expansion law if the nebula is bipolar and the slit slightly off-center. 
The third example shows different ways to produce line profiles that could be attributed to a turbulent velocity field while 
there is no turbulence in the model. 



1. THE NEED FOR 3D 
PHOTOIONIZATION CODES 

Most ionized nebulae have complex morphological 
structures. This is true not only for HII regions (Kim 
& Koo 2003, Reich et al. 1990, O'Dell 2001,Martm- 
Hernandez et al. 2003), but also for planetary nebulae 
(see many contributions to this book). High spatial 
resolution images and since a few years Integral Field 
Units observations are giving more information about 
the morphologies and kinematics of HII regions and 
PNe. 

To model such aspherical structures, a 3D photoion- 
ization codes is needed. A few of them now exist (Baess- 
gen et al. 1995, Gruenwald et al. 1997, Ercolano et al. 
2003, Wood et al. 2004). But there are some practical dif- 
ficulties with 3D model fitting of real nebulae. One is the 
enormous parameter space to consider, especially in the 
case of complex morphologies. Another one is the fact 
that 3D photoionization models are very CPU-time and 
memory consuming: several hours are needed on clusters 
or supercomputers to make a model (even if the time and 
memory can be reduced when modeling simpler geome- 
tries). 

Also, the difficulty in visualising the results does not 
help. 



2. THE TOOLS NEBU_3D, VISNEB_3D 
AND VELNEB_3D 

In this paper, we present a pseudo-3D photoionization 
code (NEBU_3D), which models aspherical planetary 
nebulae in a few minutes. The code uses various runs 
of a ID photoionization code (NEBU, see Morisset & 
Pequignot 1996, Pequignot et al. 2001), changing the 
nebular density distribution at each run according to an 
angular variation law. The set of the ID outputs (such 
as the electron temperature, the emissivities of selected 
lines at each radial position, etc.), which is defined 
according to the goal of the model, is then read and 
interpolated on a predefined coordinate cube. 

Its associated visualisation tool, VISNEB_3D, allows 
one to easily plot diagrams that are useful for understand- 
ing the physical conditions inside the models and dia- 
grams that can be directly compared to observations. The 
resulting emissivity cubes (one for each line of interest) 
can then be rotated and a projection on the sky plane is 
done to obtain surface brightness maps. 

With such a procedure, one can easily construct mod- 
els of ellipsoidal or bipolar nebulae. One can also treat 
champagne flows or even irregular nebular geometries 
provided that the photoionization source is unique and 
that the effects of non-radial ionization are not expected 
to be important. 

We also present a tool, VELNEB_3D, which can be 



applied to the results of 3D photoionization codes to gen- 
erate emission line profiles, position-velocity maps and 
3D maps in any emission line by assuming an arbitrary 
velocity field. More details can be found in Morisset & 
Stasiriska (2005). 

NEBU_3D is thus not a full 3D photoionization code 
like Mocassin (Ercolano et al. 2003). The main limita- 
tion of NEBU_3D is in the treatment of the diffuse field, 
which is not fully consistent. For example, NEBU_3D 
should not be used to model the ionization of shadows 
behind knots, as ID models cannot account in an accu- 
rate way for the ionization of the shadowed region by 
radiation produced by the surrounding material. On the 
other hand, when the radiation field is not changing dras- 
tically from one direction to the other, the use of a set of 
ID models does not lead to strong errors. The main ad- 
vantage of NEBU_3D is that it is faster by 3 to 4 orders 
of magnitude than a full 3D code when running on the 
same computer. 

More details can be found in Morisset et al. (2005) and 
Morisset & Stasihska (2005). 



3. APPLICATIONS OF NEBU_3D: A 
FEW TOY MODELS 

In the following subsections a few examples of models 
obtained by NEBU_3D are discussed to illustrate the 
capacities of 3D modelling. 



3.1. A blister and its spherical impostor 

A blister model is obtained assuming a plan parallel 
density distribution. In the case of a face-on blister, the 
surface brightness map shows a decrease of the emission 
with the angular distance to the star; the same can be 
obtained with a spherical distribution around the ionising 
star. Such a distribution is called the spherical impostor. 

The blister model presented here is a plane parallel 
slab, with density increasing exponentially with the dis- 
tance to the ionizing star, as found by Wen & O'Dell 
1995 for the Orion nebula. The gas density distribution 
law of the spherical impostor has been obtained by trial 
and error, so that the HI surface brightness is the same 
than for the blister seen pole-on. 

The Fig.^shows surface brightness maps for the blis- 
ter seen face-on (left panel), seen pole-on (center panel) 
and for the spherical impostor (right panel), for four dif- 
ferent lines. The difference between the blister and its 
impostor is perceptible for the low ionization species, 
where a ring is visible in the case of the impostor. Emis- 
sion line profiles are also different, see Morisset et al. 
(2005). 



3.2. A bipolar nebula: complexe line 
profiles 

A bipolar nebula is constructed and a linear expan- 
sion velocity law is assumed. Monochromatic images 
are shown for four lines in Fig. [2] Emission line profiles 
are obtained through various apertures and are plotted in 
Fig. [3] We can see that complex profiles (especially in 
case of narrow off-center apertures) are obtained using a 
very simple expansion law, the complexity issuing from 
the bipolar morphology. 

There are now an increasing number of instruments 
which use integral field units (IFU) to achieve spec- 
troscopy of an extended portion of the sky, either us- 
ing lens arrays, optical fibres or image slicers. These in- 
struments can thus generate 3D maps of extended ob- 
jects, with wavelength (or velocity) being the third di- 
mension. Such examples can be found in Ambrocio-Cruz 
et al. (2004) and Vasconcelos et al. (2005). Our tool is 
well-suited to also produce such maps for photoioniza- 
tion models of asymmetric nebulae. One such example 
is shown in Fig. |4j which represents the bipolar nebula 
model. 



3.3. An ellipsoidal nebula mimicking 
turbulence 

Recently, Neiner et al. (2000) and Gesicki et al. (2003) 
have claimed to find evidence for turbulence in certain 
planetary nebulae, especially in those ionized by [WC] 
type central stars. 

Two models were constructed, the first one being a 
spherical nebula with linear expansion law and high tur- 
bulence, the second one being an ellipsoidal nebula seen 
pole-on with a linear expansion law and no turbulence. 
With such a configuration, the weight of the zones of 
nearly zero radial velocity becomes large, and high ve- 
locity wings are produced as well. 

Therefore, nebulae that would correspond to such 
models would be impossible to distinguish in practise, if 
using observational setups such as those represented by 
the last two rows of Figs. [5] and [S] Note, however, that 
with a smaller slit, the two models can be distinguished, 
at least if the slit is perfectly centered (compare the sec- 
ond rows of Figs.[5]and |6j because in that case the split- 
ting of the emission lines due to crossing two emitting 
regions on the line of sight. 



4. FUTURE APPLICATIONS 

NEBU_3D can easily be adapted to any other ID pho- 
toionization code, a version using Cloudy (Ferland et al. 
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FIGURE 1. Synthetic Hj8, He I 6678A, [O II] 3726,29A and [O III] 5 007 A surface brightness maps from upper left to lower 
right. Left panel: the blister model seen tangentially (in this case the maps for [O II] 3726,29 A and [O III] 5007 A are upside down); 
central panel: the blister model seen face on; right panel: the spherical impostor model. 
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FIGURE 2. Monochromatic images obtained for three emis- 
sion lines (H/3, He II 4689A, [N II] 6583A, and [O III] 5007 A), 
for the bipolar nebula. The apertures used to compute the emis- 
sion line profiles shown in Fig. [5] are superimposed to the 
images: the sizes of the apertures are: l"x3" and 2"x2". The 
largest aperture (10"xl0") covers the entire nebula and is not 
shown. The centered apertures are drawn with solid lines. Off- 
center apertures are drawn with dashed lines. Surface bright- 
ness units are arbitrary, but the same for the four images. 



1998) is in development and a Cloudy_3D will be public 
within a few months. 



ACKNOWLEDGMENTS 

G.S. is grateful to the Instituto de Astronomia, UNAM, 
Mexico, for hospitality and financial support. The com- 
putations were carried out on a AMD-64bit computer 
financed by grant PAPIIT 1X125304 from DGAPA 
(UNAM,Mexico). CM. is partly supported by grant 
Conacyt-40095 (Mexico). 



REFERENCES 

1. Ambrocio-Cruz P., Laval A., Rosado M., Georgelin Y. P., 
Marcelin M., Comeron R, Delmotte N., Viale A., 2004, AJ, 
127, 2145 

2. Baessgen M., Diesch C, Grewing M., 1995, A&A, 297, 
828 

3. Ercolano B., Barlow M. J., Storey P. J., Liu X.-W., 2003, 
MNRAS, 340, 1136 

4. Ferland G. J., Korista K. T, Verner D. A., Ferguson J. W., 
Kingdon J. B., Verner E. M., 1998, Public, of the Astron. 
Soc. Pac, 110, 761 

5. Gorny S. K., Schwarz H. E., Corradi R. L. M., Van Winckel 
H, 1999, A&AS, 136, 145 

6. Gesicki K., Acker A., Zijlstra A. A., 2003, A&A, 400, 957 

7. Gruenwald R., Viegas S. M., Broguiere D., 1997, ApJ, 480, 
283 

8. Kim K., Koo B., 2003, ApJ, 596, 362 

9. Manchado A., Guerrero M. A., Stanghellini L., Serra-Ricart 
M., 1996, Bulletin of the American Astronomical Society, 
29, 732 

10. Martin-Hernandez N. L., van der Hulst J. M., Tielens 
A. G. G. M., 2003, A&A, 407, 957 

11. Meixner M., Kastner J. H, Balick B., Soker N., eds, 2004, 
Asymetrical Planetary Nebulae III No. 313 in ASP Conf. 
Ser. 

12. Morisset C, Pequignot D., 1996, A&A, 312, 135 

13. Morisset C, Stasihska G, 2005, MNRAS, submitted 

14. Morisset C, Stasihska G, Pena M., 2005, MNRAS, 360, 
499 

15. Neiner C, Acker A., Gesicki K., Szczerba R., 2000, A&A, 
358, 321 

16. O'Dell C. R., 2001, AJ, 122, 2662 

17. Pequignot D., Ferland G, Netzer H, Kallman T, 
Ballantyne D. R, Dumont A., Ercolano B., Harrington 
P., Kraemer S., Morisset C, Nayakshin S., Rubin R. H, 
Sutherland R., 2001, in ASP Conf. Ser. 247: Spectroscopic 
Challenges of Photoionized Plasmas Photoionization Model 
Nebulae, p. 533 

18. Reich W., Fuerst E., Reich P., Reif K., 1990, A&AS, 85, 
633 

19. Vasconcelos M. J., Cerqueira A. H, Plana H, Raga A. C, 
Morisset C, 2005, ArXiv Astrophysics e-prints/0506315 




FIGURE 3. Some physical parameters, surface brightness 
distribution and line profiles for the bipolar nebula. Upper row: 
radial distribution of the hydrogen density (left panel), the 
electron temperature (middle panel) and the expansion velocity 
(right panel). Solid curves correspond to the polar direction, 
dotted curves to a perpendicular one. Second row : from left to 
right, surface brightness distribution in the H/3, [N II] 6583A 
and [O III] 5007 A lines along the polar axis (solid curves) and 
in the perpendicular direction (dotted curves). Last three rows: 
from left to right, line profiles H/3, [N II] 6583 A and [O III] 
5007 A lines, through apertures represented in Fig. [2] (the size 
of the aperture is specified on top of each plot). Solid curves 
correspond to centered apertures, dashed curves to off-centered 
apertures. Intensity units are arbitrary, but the same for all the 
plots. 




FIGURE 4. A 3D map of the bipolar nebula model. Each 
quadrant corresponds to a different emission line. The surface 
brightness is represented by levels of grey. Superimposed on 
this image are the profiles of the same lines, integrated over 
areas of 5x5 pixels. 



FIGURE 5. Spherical nebula, with a low velocity Hubble 
flow expansion law and high turbulence. The setup of the figure 
is the same as for Fig. [5] 
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FIGURE 6. Ellipsoidal nebula, seen pole-on, with a Hubble 
flow expansion law and without turbulence. The setup of the 
figure is the same as for Fig. [3] 
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